Spectroscopy of the Kondo Problem in a Box 
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Motivated by experiments on double quantum dots, we study the problem of a single magnetic impurity 
confined in a finite metallic host. We prove an exact theorem for the ground state spin, and use analytic and 
numerical arguments to map out the spin structure of the excitation spectrum of the many-body Kondo-correlated 
state, throughout the weak to strong coupling crossover. These excitations can be probed in a simple tunneling- 
spectroscopy transport experiment; for that situation we solve rate equations for the conductance. 
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An impurity spin S coupled anti-ferromagnetically to elec- 
trons is one of the most fundamental many-body systems, the 
Kondo problem 1 1]. In the classic version, a S — 1/2 impu- 
rity is coupled with anti-ferromagnetic coupling J to an in- 
finite reservoir of conduction electrons with density of states 
p. At low temperatures, the spin is screened: its entropy is 
quenched by the formation of a singlet with the conduction 
electrons below the Kondo temperature, Tk, given approxi- 
mately by Tk ~ D e^ 1 ^ p for small values of J, with D the 
bandwidth of the conduction electrons. 

When the size of the electronic reservoir is finite, as, e.g., 
in ultra-small metallic particles 12, ufl, electron-hole excita- 
tions cost a minimum energy of order the mean level spacing 
A = 1/p. When A ~ Tk, the level spacing interferes with 
the screening of the impurity, leading to new physics |3|]. A 
finite electronic reservoir can also occur in mesoscopic semi- 
conductor structures 001: a small quantum dot acts as a spin 
5 = 1/2 impurity and is coupled to a large quantum dot with 
a finite level spacing comparable to Tk- In these latter struc- 
tures, the ratio Tk /A can be tuned by varying the voltage on 
the gates separating the two dots. 

In this paper we elucidate the manifestations of the strongly 
correlated Kondo state in the excitation spectrum of these 
finite size systems. The excitations can easily be obtained 
through measurements of the tunneling spectrum in both 
semiconductor systems j^t and metallic nano-particles 01 ■ 
We provide an exact theorem on the ground state spin of the 
system for arbitrary coupling, and make specific statements 
on the excitation spectrum. Tuning Tr-/A provides a way to 
study the parametric evolution of the system, a standard topic 
in mesoscopic physics 01 but now in a many-body context, 
through the full weak to strong coupling crossover. In the spe- 
cial case of uniform level spacing |3], we find the parametric 
evolution of the first excitation energy, SE — A F(TV/A), 
where the universal function F depends only on the parity of 
electrons in the reservoir. 

For concreteness, we consider the double-dot set-up shown 
in Fig.n in this set-up a large dot plays the role of the reser- 
voir (R) while the small dot provides the spin (S). We assume 



that the charging energies of both the large and small dots are 
large and that both are tuned in the Coulomb blockade regime. 
In addition, suppose that the level spacing of the small dot is 
large while that of the large dot, A, can be small. Under these 
conditions, the large dot acts as a finite conduction electron 
reservoir with a fixed number of electrons, N. The leads in 
Fig. \l\ are just used to probe the excitation spectrum of the 
large dot and are assumed to be very weakly coupled. 

If the tunnel coupling between the large and small dots is 
small, we can represent the small dot through a spin 5=1/2 
operator, and write down a simple model to describe the 
reservoir-spin system: 

H RS = e « c Lc Q a +JS- s(0) + E C (N - n g ) 2 . (1) 



The first term describes electron-hole excitations in the reser- 
voir (large dot), with the operators c) aa creating one body 
eigenstates on R with energy e a , wave function <\> a (r) and spin 
a. The second term describes the antiferromagnetic exchange 
interaction between the two dots, with s(0) = \f\ a ~S * a(J i fo a 
the spin density in the large dot at the tunneling position r = 
and /q ct = J2 a ^a(0) c l- F° r completeness we retain with the 




FIG. 1 : (color online) A double dot system, coupled very weakly to 
leads ("LI" and "L2"). In the Coulomb blockade regime, the small 
dot "S" behaves as a spin S that is coupled to a finite reservoir "R" 
provided by the large dot. The leads are used to measure the excita- 
tion spectrum of the system. 
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last term the electrostatic energy of the reservoir; n g and Eq 
are the dimensionless gate voltage applied to the large dot and 
its charging energy. 

Ground state spin ofH^s — We now prove an exact theorem 
for the ground state spin of the system described by Eq. Q. 
To simplify the ensuing discussion, we divide i?R_s into di- 
agonal and off-diagonal parts J/r_s = + Hon, with 

H D = JS z s z {{)) + J2 <*i fUi* + E c{N - n g ) 2 (2) 

i 

Hod = ^[S+s-(0)+h.c] -£(M flf Wa + h.c.) . 

i.a 

We first made a unitary transformation to represent the kinetic 
part of -Hr-s by a one-dimensional chain (tridiagonal matrix). 
The coefficients a*, ti and the basis states fj are fixed by the 
choice of the first basis state, /q ct . Absolute value signs may 
be used because the signs of the hopping coefficients ti and J 
can be fixed arbitrarily by choosing the phases of the f ia and 
making a tt rotation about the z axis defined by S z . 

Consider now the following set of many body states that 
span the Fock space of Eq. 0, 

\K) = (-ir^/^T-ZnT^I-^ J°) ® M > W 

with m = S z and the site labels (positive integers) ordered as 
ii < ■ ■ ■ < i _/V| and j± < ■ ■ ■ < jjy. . In this basis Hu has 
only diagonal matrix elements and Hob has only off-diagonal 
ones. Note that the Hamiltonian conserves both the total num- 
ber of electrons, N = iV-p + iVj_, and the z-component of the 
total magnetization, S z ot = (Af — N\ )/2 + m. 

It is easy to show that in this basis all off-diagonal matrix 
elements of the Hamiltonian are negative. Furthermore, re- 
peated application of H^s connects every basis state within a 
subspace (Sf ot , A) of fixed S^ ot and charge N, These facts 
allow us to generalize the proof of the "Marshall sign theo- 
rem" H to this case with several implications: the expansion 
of the ground state in the basis (0 has strictly positive coeffi- 
cients, and hence the ground state in each (S t 2 ot , A) sector is 
unique. 

These results allow us to make specific statements about the 
ground state spin of Hr_$. First consider S = 1/2. Since the 
Hamiltonian has SU(2) symmetry, we can label its eigenstates 
by S z ot and the total spin S to t • For a fixed realization of the 
reservoir R, we can always make J small enough that simple 
degenerate perturbation is convergent. With A odd, it is then 
easy to show that the ground state \G)j^o is a singlet with 
Stat = S? ot = 0. Because both \G) j^q and \G)sf ot =o,j, 
the ground state in sector Sf ot = at arbitrary J, can be 
expanded in the basis (3) with strictly positive coefficients, 
they must have non-zero overlap, implying that \G)s* =o,j 
also has total spin Stot = 0. Furthermore, this state must be 
the overall ground state, since any multiplet with Stot > 
necessarily has a state with S t z ot = 0. 

Thus we conclude that the ground state is always a spin 
singlet for any A odd. These considerations can be easily 



extended for A even to prove that the ground state spin of 
Hr.s with S = 1/2, is always Stot = for A odd and Stot = 
1/2 for A even. 

Note that so far we have made no assumptions regarding 
the energies e a and wave functions </> Q (0) and so, remarkably, 
the above theorem holds regardless of how disordered/chaotic 
the electrons in the reservoir are. We pause here to point out 
that the theorem applies to arbitrary spin, S > 1/2, or fer- 
romagnetic coupling, J < 0. It is also straightforward to 
generalize it to the Anderson model, which then takes into ac- 
count charge fluctuations on the small dot |8]. Interestingly, 
for 5=1 and J > 0, we find that the ground state can never 
be a singlet: for N odd it is a doublet while for N even it is 
a spin triplet |8]. In the 5 = 1/2 ferromagnetic Kondo prob- 
lem, we find that the ground state is always either S = 1 or 
5=1/2 depending on the parity of N. 

Excited states — In addressing the nature of the excited 
states, we now specialize to the usual S = 1/2, J > 
case, but similar arguments can be constructed for S > 1/2 
or J < 0. For an arbitrary J/ A in the S = 1/2, J > case, 
Marshall's sign theorem combined with a weak coupling ex- 
pansion only guarantees that the first excited state is either 
Stot = or 1 (1/2 or 3/2) for the odd (even) case. There- 
fore, we shall examine the excited states by perturbation the- 
ory in the two limits A > T K ( J < A) and A < T K to gain 
more insight into the structure of the excitation spectrum. 

In the weak coupling regime, A ^> Tk, given a realiza- 
tion of the reservoir R, we can always make J small enough 
that the spectrum can be constructed through lowest order de- 
generate perturbation theory. For N odd, the ground state is 
a singlet and the first excited state is a triplet with excitation 
energy 5E = K7\4> ata -p (0)| 2 A, where at op is the topmost 
singly occupied level of the large dot. The next excited state is 
a singlet of energy ~ A separated from a triplet by a splitting 
~ J . For N even, the ground state is a doublet and the first 
excited state is an 8 -fold degenerate multiplet for J = that 
gets split into two Stot = 1/2 doublets and one Stot = 3/2 
quadruplet. In general the two doublets have lower (though 
unequal) energy than the quadruplet. 

For Tk 3> A, on the other hand, the impurity spin S 
is screened by the conduction electrons, and we can use 
Nozieres' "Fermi-liquid" theory jjj. For N odd (even) we 
end up effectively with an even (odd) number of quasiparticles 
that interact with each other only at the impurity site through 

a repulsive effective interaction C/fl ~ (A 2 /Tk) ^| n l- The 
quasiparticles have the same mean level spacing A as the orig- 
inal electrons, but the spacing between two quasiparticle lev- 
els is not simply related to the spacing of the original levels in 
the chaotic quantum dot. For N odd, we find that the ground 
state is a singlet (as expected from our theorem), and the ex- 
citations start at energy ~ A. The first two excitations consist 
of a spin Stot = 1 followed by a Stot = whose energy is 
split due to the residual quasiparticle interaction by an amount 
~ A 2 /Tk. A similar analysis for the N even case shows that 
at strong coupling the ground state and first two excitations 
are all doublets. 



3 



5E 



A S= 1 



s = o 



a S - 3/2 



5E 



,''S = 1 

t 

s = o 




S= 1/2 



(N = odd) 



(N = even) 



FIG. 2: (color online) Schematic illustration of the energy eigenval- 
ues of the double dot system Eq. Q as a function of the coupling 
J for N odd (left) and even (right). In the double-dot experiment 
proposed in the text, the y-axis is like Vbias and the x-axis is like 
Vpinch- The trajectory of the excitations will show up as peaks in 
the differential conductance G. 



Remarkably, the ordering of the Stot quantum numbers of 
the ground state and two lowest excitations is the same in both 
the Tk 3> A and 7k <C A limits. It is therefore quite natural 
to assume that the order and quantum numbers are indepen- 
dent of the precise value of 7k/ A. Thus we arrive at the 
schematic illustration in Fig. [2] 

Scaling for N odd — To make further progress, we consider 
the standard simplification: |0q,(O) | 2 = 1 and e Q +i — e a = A 
independent of a [3]. Under these conditions, the spectrum, 
and thus the splitting between the ground and first excited 
states, must be a universal function of 7k /A that only de- 
pends on the parity of the number electrons on the dot. For 
N odd, for example, the singlet-triplet gap must scale as 
<5-Est/A = F(Tk/A). At strong coupling, the behavior 
F(i > 1) « 1 follows from the arguments above. To obtain 
the scaling behavior of F at weak-coupling, we sum the lead- 
ing logarithmically divergent series in the perturbative expan- 
sion of SE ST in J; this gives SE ST = J /(I - % ln(L»/A)), 
whence F(x l/ln(l/x). 

We find F(x) in the whole crossover regime through con- 
tinuous time quantum Monte Carlo calculations. The method 
we use relies on a mapping onto an XY spin chain with an im- 
purity attached to it 11(1 ll ill . The updates maintain the num- 
ber of particles (canonical ensemble). Details of the method 
and additional results shall be presented elsewhere 1 8]. 

To extract 6 Est, we measure the fraction P of states with 



S^ ot 2 = 1 visited in the Monte Carlo sampling at temperature 
T. For a fixed J and large (3 = 1 /T, P(f3) can be excellently 
fit to the form 2/(3 + e' 3 SEsT ) valid for a two-level singlet- 
triplet system (inset to Fig. [5). Repeating this procedure yields 
SEst for a variety of J and A. Plotted as a function of Tk/A 
(see Fig.[3j, we find a nice data collapse for different combi- 
nations of the bare parameters over eight orders of magnitude 
in Tk/A. These calculations thus illustrate in this special case 
that the ground state is a singlet for N odd and that the first 
excited state remains indeed a triplet for all values of J . 

Exchange coupling — So far, electron-electron interactions 
on the large dot have only been incorporated through the last 
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FIG. 3: (color online) Universal form of the singlet- triplet gap (N 
odd). Note the excellent data collapse upon plotting the extracted 
gaps as a function of Tk/A. Inset: Pqmc(/S) for D/A = 41. Cir- 
cles are data for fixed J = 0.05, 0.10, 0.13, 0.15, 0.19, 0.23, 0.26, 
and 0.30 (right to left). Lines are one parameter fits to the two state 
single-triplet form. 



term of Eq. Q. However, exchange interaction between the 
electrons also induces a Hund's rule coupling, H ex — — Jg 
(Sr is the total electronic spin in the reservoir), which will 
compete with the Kondo effect fl^l . The form of this cou- 
pling invalidates our exact theorem. However, since the cou- 
pling Js > is typically smaller than the level spacing in 
semiconductor systems 11311 and its effect (to our knowledge) 
has not been observed yet in the metallic grains we expect 
that it would only slightly modify the predicted spectrum, in- 
troducing a small splitting between the states having different 
spins in Fig. [2] There is a statistical chance that mesoscopic 
fluctuations of the level spacing will completely destroy this 
picture for some realizations of R fl4ll . 

Tunneling spectrum — In the setup of Fig.^ f° r very weak 
tunneling to the leads, the differential conductance between 
the leads as a function of bias voltage allows one to extract 
directly the excitation spectrum of the S-R system in its final 
state 1 2]. Assuming, e.g., that the reservoir initially has an 
even number of electrons, it is possible to measure the N odd 
excitation spectrum. We have predicted (Fig.|2jl that this spec- 
trum changes in an essential way as a function of the pinch-off 
voltage between the two dots that drives the system from the 
weakly coupled limit to the strongly correlated Fermi liquid 
state. A further probe is achieved by applying a parallel mag- 
netic field (Zeeman) that couples to the system via B%S* ot , 
thus splitting the spin multiplets. 

To find the conductance in such a tunneling experiment, 
we have solved the rate equations for transferring an electron 
from lead 1 to the reservoir and then to lead 2 I2l 1 1 511 - We 
assume that (1) the coupling of the lead to each state in R 
is the same (mesoscopic fluctuations are neglected), (2) the 
Kondo correlations that develop in S-R do not affect the ma- 
trix element for coupling to the leads, (3) there is a transition 
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FIG. 4: Conductance through the S-R system of Fig. Q as 
a function of the bias between the two leads for Bz — 
(0, 0.05, 0.1, 0.15, 0.2)5E (from top to bottom). For clarity, all the 
zero field splittings 5E of the resonances shown are taken equal. See 
text for more details. The initial state of the reservoir has N even 
(odd) for the left (right) panel. Both the total spin and energy split- 
ting (as displayed schematically in Fig.|2j of the device as a function 
of J can be extracted from such data. 



rate A re i that provides direct thermal relaxation between the 
eigenstates of S-R with fixed N, (4) the electrons in the lead 
are in thermal equilibrium, and (5) the temperature T is larger 
than the widths ri,T2 of the S-R eigenstates due to LI and 
L2. With these assumptions, the rate equations and their so- 
lution are standard. Results are shown in Fig.|4]in the case of 
parameters: T 2 /T 1 = 0.1, T/SE = 0.1, and A rc i/r 2 = 1. 
Note in particular two features. First, the first peak in the 
even-^odd case does not split in field while for odd^even 
it does; this provides a definite signature of the parity of N. 
Second, for even^odd, the behavior of the excited states as a 
function of field is quite different: the triplet excitation splits 
into two symmetric peaks while the singlet excitation splits 
very asymmetrically. Thus the singlet and triplet excitations 
can be identified unambiguously and so the behavior shown 
in Fig.|2]verified. 

Larger spin — Our theorem and the analysis of excited states 
can be easily generalized to the underscreened S > 1/2 sin- 
gle channel Kondo problem, which has indeed been realized 
experimentally 11611 . Interestingly, for N even and S = 1, it 
follows from our theorem that the ground state for all Xk/A 
has S*tot = 1. At weak coupling, we deduce this simply from 
perturbation theory and also that the first excited state is a sin- 
glet at energy ~ A. In the opposite limit Tk ^ A, the flow 
is to the usual underscreened fixed point jjj: one is left with a 
5 = 1/2 impurity coupled/e rromagnetically to an odd num- 
ber of quasiparticles in the reservoir. Since the ferromagnetic 
Kondo problem flows naturally to weak coupling we are 



again justified in doing perturbation theory, and we recover 
that the ground state Stot — 1 is separated from the first sin- 
glet excitation by an asymptotically small energy. Thus, the 
signature of the underscreened Kondo effect is a reduction of 
the gap between triplet ground and singlet excited states, of 
logarithmic form as a function of Tk / A, as one tunes from 
weak to strong coupling. 

In conclusion, with a combination of an exact theorem, per- 
turbative arguments, and Monte Carlo simulations, we have 
made specific predictions on the low-lying excitation spec- 
trum for a double dot system and a metallic grain contain- 
ing a single magnetic impurity. The excitation spectrum was 
found to reflect the formation of the Kondo-correlated many- 
body state in an essential way. Further by numerically solving 
the corresponding rate equations, we have shown that the gate 
voltage and magnetic field dependencies of conductance mea- 
surements can be used to observe our predictions directly in 
an experiment. 
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